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Quantum Field Theory is applied to study an electron plasma under an intense 
neutrino flux. The dispersion relation of the longitudinal waves is derived and the 
damping rate is calculated. It is shown that in the case of Supernova emission the 
neutrinos are not collimated enough to cause plasma instabilities associated to a 
strong neutrino resonance effect. 



1 Introduction 

Bingham et. al. have studied the behavior of an electron plasma under a 
neutrino flux and concluded that neutrino fluxes as intense as the ones pro- 
duced in Supernovae cause plasma instabilities. If true this could provide the 
physical mechanism of energy transfer from the neutrinos to the medium that 
would explain the Supernovae explosion. The weak interaction between the 
neutrinos and the plasma was described with the concept of a ponderomotrive 
force acting on the electrons, 

F pond = -V2G F d v Vn v , (1) 

proportional to the gradient of the neutrino number density (c' v = ±1/2 + 
2sin 2 #W; (+) f° r v e arL d (— ) for z^, v T ). The neutrino wave function was 
assumed to obey a naive Klein-Gordon equation modified to include matter 
effects. However, that equation does not account for the neutrino-spin degrees 
of freedom and those works were subject to justified controversyEm. 

More recently, a classic .kinetic theory was applied to describe the elec- 
tron and neutrino dynamics □. In that work the neutrinos too are subject to 
a ponderomotrive force analogous to Eq. (Q) with n v replaced by the electron 
number density n e . Again the lepton spin and chiral structure of the weak 
interactions were ignored. In addition, the dispersion relation derived for the 
plasma waves is substantially different from the one previously obtained with 
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the Klein-Gordon equation El. Nevertheless, the authors reiterated the claim 
that the neutrinos produce resonant effects for certain modes of plasma oscilla- 
tion which cause instabilities characterized by large growth rates proportional 
not to G F n v but to a smaller power of that number. 

The aim of the present work was to apply consistent Quantum Field The- 
ory techniques to obtain a kinetic description consistent with the quantum 
mechanics and spin content of particles and interactions. First, we derive 
the dispersion relation of the longitudinal waves of an electron plasma in the 
presence of a neutrino flux and afterwards, analyze the possibility of neutrino 
induced plasma instabilities. 



2 Classic Kinetic Theory 

It proves useful to go first to a classic description before going to Quantum 
Field Theory. The system of interest is a large number of electrons, neutrinos 
and anti-neutrinos. Classically they are described by distribution functions, 
f e , /„ and fp, on the respective position and momentum phase spaces. In low 
dense plasmas the collisions are less important than the collective interactions 
and the evolution of the distribution functions f(t, x,p) reduces to the Vlasov 
equations Era, 

f < 2 > 

The dynamics is defined by specifying the velocity v and force F as functions of 
(i, x,p) for each of the particles species. In the electron case v=p/ (p^+ml) 1 / 2 , 
but the neutrino mass will be neglected. 

At low energies the electron interactions are described by the Standard 
Model effective LagrangianLl 

Ant =e^eYe- V2G F (P L -/^ L ) e^(c' v ~ c' A -/ 5 )e , (3) 

where e is the positron charge, G F the Fermi constant, and c' v — c' A + 2 sin 2 #w 
with c' A = +1/2 for v e and c' A — —1/2 for u^, v r . By analogy with electro- 
dynamics one immediately infers the Lagrangian of a classic spinless electron 
as 

L e =(eA^-V2G F c' v J^)x^. (4) 



The average current density J£ — (^l7 m ^l), equal to the difference between 
the neutrino and anti-neutrino current densities, Jjf = j£ — j£, acts in the 
same way as the electromagnetic potential A^. Hence, the electro- weak force 
applied to an electron is a straightforward generalization of the Lorentz force, 

F e = -e(E + v e AB) + V2G F c' v (E we + v e A B we ) , (5) 
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with weak-electric and weak-magnetic fields given by 



dJ - 

E we = -VJ° - , B we = V A J v . (6) 

In a similar fashion, the weak forces applied to massless neutrinos is 

F v = V2 G F {E wv + v u l\ B m ) (7) 
for neutrinos (vl) and —F v for anti-neutrinos (vr), where 

dJ^e - - 

E WI/ — — V J we -^j— , B wu = V A J we , (8) 

are obtained from the weak electron current, 

JL = (ej^ic'v - c'a 75) e) = c' y J e M - c' A Jg e . (9) 

It contains an axial current as well, that is different from zero if the electrons are 
polarized: for an electron density n e and average polarization (<?), (^775 e) = 
(a) n e . 

The first point to noteJ.s that the weak forces contain other terms than 
the ones identified before an namely, those with the vector current densities 
J„ and J we . That is true even in the limit of non-relativistic and unpolar- 
ized electrons because the neutrinos are ultra-relativistic. Then, the following 
approximations apply: 

F e S ~eE~V2G F c' v | Vn„ - Vn p + ^ - ^ j , (10) 



Fv = -Pp = -V2G v c' v Vn e + ^f . (11) 




The electron current can be neglected, as presented in the conference, except 
for extreme wavelenghts, comparable to the inverse plasma frequency times 
the velocity of light. For longitudinal waves V A j e is still zero and does not 
contribute to the neutrino force. 

The proper modes of a system are usually investigated by expanding the 
distribution functions around certain zero order functions f°. If the back- 
ground is static and uniform, /° only depend on the momenta and 

f(t,x,p) = f(p) + Sf(t,x,p) . (12) 
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The Vlasov equations (^|) are linearized by neglecting the quadratic terms in 
Sf and the Fourier analysis is applied afterwards. Whenever the forces vanish 
at zero order, as in the case of the electro-weak forces, the Fourier transforms 
obey the equations 

df° 

(w - k-v)Sf(k,p) +iF(k,p) ■ = , 

where u> = k° is the plasma frequency and k the wave vector. The problem of 
deriving the dispersion relation is facilitated by the fact that the forces (fiol), 
( |i"T| ) depend linearly on the current densities 



f{t,x)= d*pf(t,x,p)vx (13) 

(v° = 1). Then, the fluctuations 

5r(k)= I > p ^L-^.iF(fc,pK(p) (14) 
J Lu-k-v OP 

form a closed system of linear equations whose secular equation determines the 
dispersion relation. 

The Langmuir waves U'u are electron density waves associated with an 
electrostatic field {B = 0, rot E = 0) whose potentials are 

A°(k) = ^6n e (k) , A(k)=6. (15) 
k^ 

They are also called longitudinal electromagnetic waves and sometimes plas- 
mons in the Particle Physics literature When the electro- weak forces (^p|), 
( |TT| ) are substituted in (|l4|), one obtains 5n e as a linear function of 8n e , 
8n v — drip, and 5j v — Sj p and, in turn, 5n v , 6j v , Sn p , 5jn proportional to 5n e 
(Sj e = uikSn e /k 2 ). It yields an equation of the type Sn e = —xSn e , where the 
susceptibility \ nas > m addition to the purely electromagnetic x E m > a w eak 
interaction term \ w ■ The plasma dispersion relation is 

l + X BM (u,k)+X w {",k) = . (16) 

Xem on ly depends on the electron charge and distribution function. It is shown 
in many books that (e 2 = Ait a) 

X EM = i ffp—^yji-k. (17) 
k 2 J oj- k-v e op 
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In the presence of a neutrino beam, the weak correction is given by 



,2X EM {1 u 2 \ [«k-(k-uv u ) (dfl t <9/<T 



Xw . - 1Wf {1 - ¥) J fp-Z^ [if + (18) 

Notice that despite the relative sign in the expression of the electron force, 
Eq. (|l^), the v and v beams contribute in the same way to the plasma sus- 
ceptibility. The reason is the v and D density fluctuations get also opposite 
signs from the forces F v and F„. That expression of \ w differs from the result 
obtained inta by the term —tx><o v under the integral and the factor 1 — Lu 2 /k 2 , 



which come from the time derivative of J„ and j e in the forces ( pL0| ) and (11) 
respectively. 



3 Quantum Field Theory 

A basic property of the weak interactions is their chiral structure i.e., the 
differentiation between the left and right-handed components of the fermions. 
This aspect is ignored in the classic theory because it does not include the 
dynamics of the spin degrees of freedom. At the fundamental level the basic 
entities are the quantum fields but in the situation of interest they interact 
within a bath of plasma and neutrinos. For a medium in thermal equilittj 
rium there are well established Finite Temperature Field Theory techniques □ 
however, the passing neutrinos are not in thermal equilibrium and the plasma 
waves represent themselves deviations from equilibrium. The matter content 
is represented by a density of states p. In the Heisenberg picture, p is constant 
and the fields evolve in time. In the Standard Model they are subject to local, 
relatively simple interactions and we caa. take advantage of that to solve the 
field equations of motion perturbativelyH. 

In matter, the one-particle propagators are different from the propagators 
in vacuum. I define the distribution functions as a difference between the 
expectation values in matter and in vacuum of the time-order products of the 
fields in the Heisenberg picture i.e., 

(Tr(x)^(y)) p = (o\Tr(x)^(y)\o)-r p (x,y) ■ (19) 

That f a P(x,y) represent one-particle distribution functions can been seen by 
evaluating observables like the current density 

J»(x) = {i>{x)^ip{x)) - (o|$(a07"t/>(a:)|Cl) (20) 
= tr{f(x,x)^} , 
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where the trace is over the spinor indices. 

In absence of interactions one knows the exact expressions of the quantum 
fields and exact distribution functions can be easily specified. The function 

fe(x, = jjil /e(P) 0>e + ™e) e^'^ (21) 

describes a homogeneous (invariant under translations) distribution of unpo- 
larized electrons with momentum distribution f®{p), as shows the expression 
of the current density, 

jf{x) = tr {f e (x, x) r} = J d 3 pf e (p) J . 

Likewise, uniform beams of massless neutrinos have a distribution function 
(associated to the field vl) 

f v (x,y) = ^l J^L (f^fre-^-rt - f^fre+^-y^ , (22) 

where /° and f§ represent the (left-handed) neutrino and (right-handed) anti- 
neutrino momentum distributions. f®{x,y) and f®{x,y) will be used as back- 
ground distribution functions which implies a plasma with no positrons and 
therefore temperatures below the electron mass. 

The functions f~ (x, y) and f®{x, y) stop to be consistent with the definition 
( |l9| ) when the field equations of motion include interactions. The deviations 

5f(x,y) = f(x,y)~f°(x,y) (23) 

can be calculated perturbatively i as part of higher order corrections to the 
one-particle Green functions ( |l9| ) along with other radiative corrections such 
as vacuum self-energies and matter induced 'potential' energies. Sf e and 5f v 



fe fe 




Figure 1: Diagrams contributing to <5/„. 

G 




Figure 2: Diagrams contributing to 8f e . 

are the corrections associated to the very matter fluctuations 8f(x,y) and, in 
first approximation, they are given by the diagrams in Figs. 1 and 2. 

At one-loop, the Fourier transforms of 8f(x,x), Sf e (k) and Sf u (k), are 
just linear combinations of the leptonic current fluctuations namely, 8J£ — 
tT{8f v ^}, 5 J? = tr{<5/ e7 ^} and SJg e = tr{8f e ^ 5 }. As a result, 6Jg(k), 
5J^ e (k) and 8J£(k) form a closed system of coupled equations, 8 J — AS J in 
a matrix notation. In principle one could derive the dispersion relation u>{k) 
from the condition det(l — A) ~ but that is not easy. In particular, the 
Lorentz structure is non-trivial because the plasma rest frame, wave vector 
k and neutrino fluxes define privileged directions in space-time. The first 
possible simplification is to take the macroscopic limit by neglecting u> and k 
with respect to the one-particle energies E e , E v , making in particular (p e ± 
k) 2 — ml pa ±2k ■ p e and (p„ ± k) 2 sa ±2k ■ p u in the electron and neutrino 
propagators. It turns out that in this limit the field theory results for the 
currents 8J£ and 8J£ are equivalent to the classic expressions ( |l4| ) of 8jg 
and 8j£ — 8j£ respectively, with forces F e (k,p) and F u (k,p) given exactly by 
the Fourier transforms of the forces (Q) and (0), derived from the appropriate 
classic theory. 

Of course, the axial current 8 J^ e can only be calculated within field theory. 
One expects from the Standard Model Lagrangian (0) that the neutrino vector 
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current polarizes the electron spin like a magnetic field does on a magnetic 
moment. The density SJ® e is suppressed for non-relativistic electrons but the 
axial current S Js e gets polarized along the component of the neutrino velocity 
orthogonal to the wavevector k. It is also accompained by a transverse electric 
current both of them proportional to G|. However, their contribution to the 
dispersion relation only comes at the Gp order and can so be discarded. The 
classic dispersion relation (|l8|) is recovered for the non-relativistic plasma. 

4 Plasma oscillation. Damping. 

Let Wpi(fc) denote the plasma frequency as a function of k in_the absence of 
neutrinos. For wavelengths much larger than the Debye radius Era, the electro- 
magnetic susceptibility is given by X EM (w, k) — — u^fcj/w 2 . After applying 
this relation in the Eqs. (|l^) and ( f[8| ) the dispersion relation of the non- 
relativistic plasma in the presence of neutrinos reads as 

" Upl{k) ~ 2na P {k Uj) J dp w -k-v v [dp + dp) k > 

(24) 

or, after integrating by parts, 



- _ _ „ r& (/ . + m p>0t . (25) 

P 27TQ! k 2 J E v (<jJ-k-V v ) 2 



Again, as a^consequence of the discrepancies in the weak forces, these results 
differ fromcl namely, by the factor (1 — w 2 /fc 2 ) 2 in the second equation. 

The constant c'y is much smaller for v ^ , v T (c' v ~ —0.04) than for v e (c' v ~ 
0.96), so only v e and v e deserve to be taken in consideration. The expressions 
above provide a good perspective of the real impact of the neutrinos. Keeping 
only the main factors, 

Gf n v 



2 2 2 '"v f2 

w -uj pl cx uj p1 — - — G F fc 



clearly indicates that the neutrino effect is severely suppressed by Gp. The 
only potential exception are the oscillation modes that catch the neutrinos in 
the resonance u> = k-v v . 

If the neutrino spectrum goes from one side of tha-pole to the other the 
situation is analogous to the electron Landau dampingtffl. The integral in Eq. 
( p4j ) separates into a negligible principal part and a complex quantity that is 
determined by the so-called Landau prescription, here u> — k-v v + ie, s — ► 
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+ (identifies with the Feynman boundary conditions on the propagators). 
Denoting the imaginary part of u) as — ij, the neutrino contribution to the 
decay rate 7 is then 



Hardy and Melrose a obtained this result from the calculation of the decay 
rate of neutrinos into longitudinal photons. However, by the very nature of 
such calculation the full neutrino contribution to the dispersion relation is not 
derived and second, cannot be utilized to investigate possible hydrodynamic 
instabilities. 

7w S oes as ^1 an d is exceedingly small. If however all the neutrinos were 
on the top x£ the resonance they could generate an hydrodynamic instability. 
The claim EJu was that this occurs in the conditions of Supernova neutrino 
emission causing much larger growth rates. The question I rise is, in the end 
one has to check whether or not the entire neutrino flux lies in the resonance 
i.e., whether \u> — k-v v \ stays within the calculated resonance width | | . The 
neutrino velocities only spread in direction. At radii, r, much larger than the 
neutrinosphere radius, i?„, they move essentially in the radial direction yet, the 
angular spread V ss 2R u /r is finite. Under the resonance condition, k-v v w io p \, 
the interval of variation of k-v v = \k\ cos 9 u k lies in the range 



The left-hand value holds for a vector fcwith radial direction whereas the right- 
hand bound is necessary to exclude the situation of Landau damping and the 
result @. 

The largest values of I7I are obtained by assuming that k • v v is approxi- 
mately constant over the neutrino spectrum and can be factorized out of the 
integral in Eq. ( p5| ) . One obtains growth rates of the order of 





(26) 



(27) 





for the resonant modes (hv v ~ ^pi)- 7opt. varies with tan 4 / 3 9 v k and Afc-tv with 
tan 6 v k , so the most favorable case for a hydrodynamic instability is a vector k 
orthogonal to the radial direction i.e., Ak-v^ ~ lu p \ and tan#„fc w 1/0 V . Then, 
the condition that the neutrinos be inside the resonance requires 
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Knowing how small the energies Gf n v andJ^F Ti e are, that looks quite unlikely. 
For numerical estimates I used the valuesBl L Uc — 10 53 ergs/s, E v = lOMeV, 
n e = 10 33 cm~ 3 and 9 V — lOOkm/r. For 6 V — 0.1, 7 is of the order of 10~ 7 , 
quite small indeed. Furthermore, since the neutrino density varies as 1/V 2 , 
7opt./w p i is proportional to (ngr 2 ) 1 / 3 and necessarily drops to zero for large 
enough radius. One may also fix the value of the electron density. It turns out 
that the parameter 7 approaches the unity only at a radius equal to 10 13 km 
for n e as large as 10 33 cm -3 ! That is clearly absurd. 

It should be also taken in consideration that for wavelengths smaller than 
the Debye distance (k > k^), the plasma waves are electromagnetically Lan- 
dau dampedtm That puts a further upper limit on the interesting values of 
tan 6 v k for the resonant neutrinos (< kjj / uj p i) in disfavour of the hydrody- 
namic instabilities. 



5 Conclusions 

We have shown that Quantum Field Theory yields results that are compatible 
with the classic kinetic theory provided that one uses the right specification for 
the forces. Simple considerations show that the weak forces between electrons 
and neutrinos are of the same type as the Lorentz electromagnetic force. In 
particular, the force done by a neutrino flux on a non-relativistic electron is 

F e = -V2 G F c' v (vn u + ^ . (28) 

In Quantum Field Theory the spin degrees of freedom are also taken in account. 
Under a neutrino flux the electrons get polarized at Gp order but the final effect 
on the plasma dispersion relation is of higher order and can be neglected. 

In what concerns possible plasma instabilities caused by the neutrino flux 
in a Supernova envelope star we have shown that the neutrinos are never 
collimated enough to produce strong resonance effects. As a consequence, the 
neutrino induced growth rates are of the Landau damping type, suppressed by 
Gp. 
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